Non-Galilean response of Rashba coupled Fermi gases 
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We consider tiie effect of a momentum kick on tlie ground state of a non-interacting two- 
dimensional Fermi gas subject to Rashba spin-orbit coupling. Although the total momentum is 
a constant of motion, the gas does not obey the rules of Galilean relativity. Upon imprinting a small 
overall velocity to the non-interacting gas, we find that the Fermi sea is deformed in a non-trivial 
way. We also consider a weakly repulsive Fermi gas, and find, from its Hartree shift, that the to- 
tal ground state of the system may change into a deformed, finite momentum ground state as the 
repulsion is increased beyond a critical value, without the need of any external Zeeman fields. We 
also discuss possible experimental signatures of these effects. 
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Recent progress on synthetic spin-orbit coupling for 
neutral atoms provides new ways to study its ef- 

fects in both bosonic [^-0] and fermionic d, Q systems. 
For atomic Fermi gases, many theoretical studies have 
shown that the spin-orbit coupling has nontrivial conse- 
quences on the BCS-BEC crossover when an attractive 
interaction is present Even when the attraction 

is small, the BCS-BEC crossover can be induced by in- 
creasing the spin-orbit coupling strength. Due to the 
absence of Galilean invariance in the spin-orbit coupled 
gases, fermionic superfluids can exhibit highly non-trivial 
properties These systems can undergo quantum 

phase transitions to some topological superfluid phases 
[18l421|. and it has also been shown that the finite mo- 
mentum pairing state or the so-called FFLO state can 
be energetically favored in spin-orbit coupled Fermi gases 
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26[. However, these more exotic phenomena require 
an additional, carefully engineered Zeeman field. 

In this Letter, we study the response of a Fermi gas 
with pure Rashba spin-orbit coupling to a small momen- 
tum kick, in two spatial dimensions for concreteness. We 
find that, even in the low-momentum regime, the Fermi 
sea is deformed in a non-trivial way, is highly degener- 
ate, and that its energy is bound from above by that 
of trivial Galilean-like boosts. We then consider repul- 
sive interactions. We find, using mean-field theory, that 
the ground state acquires a finite momentum for interac- 
tion strengths beyond a critical value. The Fermi sea is 
consequently deformed, without the need to include any 
external Zeeman fields. 

The single-particle Hamiltonian of the system is given 
by [27| 



p2 + fi2A2 f,x 

Ho = 1 H cr • p, 

Zm m 



(1) 



where a — [a^, <^y) is the vector of spin-1/2 Pauli matri- 
ces. The constants of motion are the helicity H = cr • p, 
with its corresponding eigenvectors IV'^^') forming the 



helicity basis, with eigenvalues /i = ±|p| for momentum 
p = ftk, and the momentum itself. The common cigen- 
states of the Hamiltonian, helicity and momentum have 
the form 



ikr 



(2) 



where is the polar angle of k, given by tani/) = ky/k^. 
The two energy branches of the system, corresponding 
to negative (lower branch) and positive (upper branch) 
helicity, respectively, are given by 



£±(k) 



2m 



(|k|±A)= 



(3) 



For the sake of simplicity, we will assume that the 
Fermi gas is dilute enough that in the ground state parti- 
cles only occupy negative helicity states. For this to hold, 
the Fermi energy must be Ep = h^kp/2m < h^X^/2m, 
with the Fermi momentum kp = irp/X. Defining the 
only dimensionless parameter of the gas z = np/X^, the 
single-branch condition reads z < 1. The ground state 
energy density £ is easily calculated and reads 



6m 



TTZp. 



(4) 



The Fermi sea in two dimensions has the form of a 
concentric annulus. For convenience, we define "radii" 
Rj ^ X — kp = A(l — z) and Rq ~ X + kp ~ A(l + z) for 
the inner and outer circumferences, respectively, which 
are the borders of the Fermi sea. 

We study now the properties of the lowest energy Fermi 
sea when the gas is given an infinitesimally small momen- 
tum kick. The momentum kick per particle is denoted as 
ko. Without loss of generality, we choose the momentum 
kick in the positive x-direction, kg = k^kx, with fco > 0. 

We begin by considering the obvious Galilean-like 
boost, consisting of assigning an extra momentum fco to 
each fcrmion in the Fermi sea. We note that, formally. 
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these are not Galilean boosts since they also involve a 
non-trivial spin rotation, as is observed from Eq. ([2|)- 
Geometrically, these transformations correspond to dis- 
placing the Fermi sea to the right by an amount . The 
excess energy density A£c is readily calculated to O(fco) 
and is given by 



(5) 



Below, we will see that the energy of these transforma- 
tions is an upper bound for the ground state energy at 
finite momentum, and is only attained for a saturated 
lower branch (z = 1). 

If the momentum kick per particle ko to the gas is 
small, the ground state in the thermodynamic limit can 
be modelled by infinitesimal transformations, which de- 
scribe the response of the system to small momentum 
kicks. These must fulfill two conditions: (i) the density 
p is preserved and (ii) the momentum per particle of the 
resulting Fermi sea equals ko. If we consider infinitesi- 
mally small momenta, we can safely rule out breaking the 
Fermi sea into disjoint pieces. We are then left with three 
possibilities, namely scalings of the radii Rq and Ri by 
a factor a and /3, respectively, displacements of the inner 
and outer circumferences, and multipolar deformations 
of these. 

Out of the three types, at small momenta scalings are 
ruled out since they do not contribute to adding mo- 
mentum but they do increase the energy. Scalings that 
preserve the density of the system are restricted to the 
following relation 



(6) 



We now consider displacements of the inner and outer 
circumferences with respect to each other, which always 
preserve the density of the system. The resulting momen- 
tum density is easy to calculate if we regard the empty 
region |k| < i?/ as a virtual Fermi sea of holes with den- 
sity Pi, and the partially filled region |k| < Rq as a vir- 
tual Fermi sea of particles with density po- If we displace 
the inner and outer circumferences by momenta q/ and 
qo , respectively, the momentum density of the system is 
given by 




koP = qopo - mpi- 



(7) 



Finally, deformations of the Fermi sea can be 
parametrized in polar coordinates by allowing the radii 
become angle-dependent as Ri{4>) = Ri + fi(4>) and 
Roi4') = Ro + fo{4>), with // and fo real periodic func- 
tions. Clearly, not all deformations preserve the density. 
This is given by p — po ~ pi, which implies the following 
condition 



dcl>{[fo{<P)f^[fi{<P)f} 



= 0, 



(8) 



or, cquivalently ||/o|| = ||//||) with || ■ || the norm on 
L^([0,27r)). This condition is very relaxed, since it al- 
lows very different deformations of the two circumfer- 
ences provided that their norms are equal. Since, by 
convention, we have chosen to imprint a momentum to 
the system in the positive x-direction, we can restrict 
ourselves to deformations in the horizontal axis. There- 
fore, we parametrize the functions fi {i = 0,I) by a 
multipolar expansion of the form 



/TT 



^ ^{t) cos(m0) 

m>l 



(9) 



where every coefficient Cm is real. The condition on the 



norms takes the form X]m>i(c^'')^ = J2 



m> 



l(Crn')^. 



The 



momentum density of the deformed Fermi sea is given by 
kop = kopkx, with 



koP 



1 



(2^)2 



27r rBo+foW 

d(j) I dkk'^ COS (f>. (10) 

JRi+fiW 



For an infinitesimally small momentum kick, we can re- 
strict the expansion in Eq. ([9|) to lowest (dipolar) or- 



der, in which case condition ^ reduces to c^^^ = ±c{ 



„(o) 



JO) 



TT, with the upper (lower) sign cor- 



Defining 

responding to equal (opposite) dipolar coefficients, we 
have fco = qt , while for c^' 



in 



SO) 



kop 



R] 



we have 



(11) 



To lowest order in ko, is given by q^ = 2zfco/(l + z^). 
Intuitively, dipolar coefficients with opposite sign will 
yield a lower energy, since a weaker deformation (|gj~| < 
\q^\) is needed in order to have the required momentum 
/cq. The excess energy density A£± for dipolar deforma- 
tions is given, to leading order in q^ , by 



±\2 



A£± 



4m 



(12) 



From the above relation, we see that for equal deforma- 
tions the excess energy is that of Galilean-like boosts, as 
in Eq. ([5]), A£+ = ASq- For opposite dipolar coeflBcients 
we obtain a lower energy given by 



AS. 



h^k^ 



l + z2 



P, 



(13) 



which is bound from above by A£q. The excess en- 
ergy due to pure deformations can be lowered further by 
considering either displacements or both displacements 
and deformations. We will show next that the result- 
ing ground-state at small non-zero momentum is highly 
degenerate. 

We now study the energy of the Fermi gas when we 
impose displacements q/ and qo of the inner and outer 
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circumferences, together with deformations of opposite 
coefRcients, parametrized by qi = q^. The energy den- 
sity of the system is calculated as the difference between 
the energy of the virtual Fermi sea of particles, So, and 
that of the virtual Fermi sea of holes, £j, where 



1 



(27r)2 



2ir 



dfcfce_(k + q,), (14) 



with i — 0,I. To lowest order in the momenta, we can 
isolate the excess energy from the energy at zero momen- 
tum, which takes the form 

(27r)2|jA£- = nkliql - qj) + TrXkpiql + qj) 

+ 2nkF [kpqiiqo + qi) + Mi{qo - qi)] 



2'KXkFqi ■ 



(15) 



If the momentum per particle is set to kg = k^kr^^ the 
three momenta of the system must fulfill the condition 



qopo - qiPi 



1 + 
2z 



pqi = kop. 



(16) 



Using the momentum constrain above to clear qj in fa- 
vor of qo and qi, and inserting this in to Eq. ([15]), we 
obtain an expression for the excess energy density, which 
is minimized for any qo satisfying 



qo 



-qi 



2z{l + z) ^ 
1 -h3z2 



(17) 



with the only constrain that both qo and qi have to be 
small (^ kp) for our theory to be vahd. We thus see 
that the ground-state at fixed finite momenta is infinitely 
degenerate. Finally, the excess energy density is given by 



m 1 + 3z' 



:P, 



(18) 



which is bound from above by ASo for all z and by A£_ 
for each z. 

We consider now the effect of weak interactions, which 
we model with a zero-range potential. At mean-field 
level, the Dirac delta needs not be regularized, and has 
the form (7(5^^^(r), where r is the interparticle distance. 
The dimcnsionlcss interaction parameter ^ = mg/ATrh^ 
[28j must here be compared with the only dimensionless 
parameter z of the non-interacting gas. The typical inter- 
action energy [i^ should be smaller than the Fermi 
energy Ep = Trh^pz/2m, which gives f_< z/2. 

The Hartree shift f^^-* is given by [29[ 
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4(27r)4 



dk / rfk' cos( 



0k')- (19) 



The first term is independent of the shape of the Fermi 
sea, and the second term vanishes identically in the non- 
interacting ground state. However, the second term re- 
duces the interaction energy for repulsive interactions 



{g > 0) if the Fermi sea is modified to acquire a non-zero 
total momentum. Without loss of generality, we again 
consider a momentum per particle kg in the i-direction, 
in which case the extra contribution to the Hartree shift 
reads 



X: 
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4(27r)^^j. 



(27r)2 



dk / dk' cos(0k 
p 

1 

dk cos ( 



') 



(20) 



In the following, we minimize the total energy excess 
A£t — AE-\-X at a constant, infinitcsimaly small momen- 
tum ko, with respect to the displacement parameter qo. 
Due to the large degeneracy of the system, we consider 
displacements only without deformations. The Fermi sea 
is the set F of momenta k defined as F = {kj^/ < |k| < 
Ko}, with = + qj + 2\qi cos (j) {i = /, O). The 

extra Hartree shift of Eq. <\2Q\ is easily calculated and 
reads 



9P 
GAttz 



[(1 + z)qo - (1 - z)qi] 



(21) 



where qo and qi are related via Eq. with q^ = 0. 

The energy excess A£t is minimized at the point 



qo 



and its value is 



^iz-l) + z + l 
^z{z^ - 1) + 1 + 3z^ 



2koz. 



A£t = -F(e, z) 



h^kl 



where 



J-(^,z) = z 



^Z(z2-l) + 1 + 3^2 



(22) 



(23) 



(24) 



As we infer from Eq. (P^ . the system's ground state, to 
first order, will change from the non-interacting Fermi sea 
to a Fermi sea with an infinitesimally small momentum 
- which denotes continuity - at a critical value ^ = £,c{z) 
where J^(^c, z) = 0. This is given by 



ic{z) = Z. 



(25) 



It is important to note that, while we are considering the 
system's non-interacting energy to order O(fco), the in- 
teraction energy is not approximate but is indeed propor- 
tional to /cq , and therefore the above critical value is exact 
at the mean-field level, and is momentum-independent. 
The critical value we have obtained is beyond the mean- 
field regime. This is a typical situation that also occurs, 
for instance, in the study of the ferromagnetic transition 
[30t in repulsive Fermi gases jsH-Hl], and is resolved by 
going beyond first-order perturbation theory, which can 



result in an apparent first-order transition |33l . |34| , or us- 
ing non-perturbative methods (ssl . 36 1. which predict a 
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second-order phase transition and are in good agreement 
with Monte-Carlo simulations 



37| 



Our findings are the 
natural starting point for higher-order corrections [2^ 
and non-perturbative treatments. Moreover, if as our re- 
sults suggest, the transition to finite momentum is con- 
tinuous, the critical point can be obtained by calculating 
only the second-order response function z) in im- 
proved treatments of the interactions. 

The low-momentum theory we have presented only 
predicts whether the system evolves towards a finite mo- 
mentum ground state in a continuous manner. This is to 
say that our theory can describe, correctly, only deriva- 
tives in the energy at zero momenta 

The first derivative of the energy with respect to fcg 
vanishes always at fcg — 0, while its second derivative 
z) > for ^ < and z) < for ^ > (strictly 
speaking, until its pole at ^oo = (1 + 32:^)/(z — z^) » ^c, 
which is far beyond the limit of validity of the mean-field 
theory). The change of sign in z) denotes the tran- 
sition from a minimum to a maximum at feo =0, which 
implies that for repulsions stronger than the critical value 

the system's ground state has non-zero momentum. To 
calculate the actual momentum of the ground state for 
^ > it is necessary to go beyond leading order in fco- 

So far, we have considered the momentum kick of the 
gas to be given in a particular direction. In an exper- 
iment, however, no direction is in principle preferred. 
Obviously, after the critical point is reached, the ground- 
state is highly degenerate, since the same momentum in 
any direction yields the same energy. The many-body 
wave function is therefore an arbitrary superposition of 
states with momenta pointing at different directions. For 
instance, we may expect to observe an equal superposi- 
tion of these states, up to arbitrary phases, which yields 
a circularly symmetric momentum distribution. Momen- 
tum distributions are especially relevant to cold atom ex- 
periments where these can be obtained via timc-of-flight 
measurements [sst . 

We can easily map out, starting from a Fermi sea 
at finite momentum in a particular direction, the afore- 
mentioned circularly symmetric momentum distribution. 
The momentum distribution is denoted by 7^(k;), and is 
such that p = J (ikn(k)/(27r)^. We define the integrated, 
angle- independent momentum distribution ri(fc) as 



2TTn{k) 



27r 



l(k). 



(27) 



Geometrically, h{k) is the arc length of occupied states on 
a circumference of radius k in momentum space. Closed 
circumferences map into unit length - fully occupied 
states - while open arcs give shorter lengths - smaller 
average occupations. The circularly symmetric momen- 
tum distribution arising from the superposition of the 




FIG. 1: (Color online) Integrated momentum distribution 
(black solid line) at critical interaction strength ^ — z, with 
fco/A = 5 • 10~^ and z = 1/5, compared to the non-interacting 
momentum distribution (red dashed line). 



different states is obtained as a surface of revolution by 
rotating the Fermi sea at finite momentum in a particular 
direction with respect to the origin in momentum space. 
Clearly, the resulting momentum distribution coincides 
with h{k). 

In Fig. [TJ we show the integrated momentum dis- 
tribution for a particular case of an interacting Fermi 
sea at non-zero total momentum compared to the non- 
interacting momentum distribution. There, we observe 
a shortening of the unit occupation plateau, together 
with a smoothcning of the momentum distribution at the 
edges of the Fermi sea, with an obvious change in con- 
cavity, which can be a relevant experimental signature 
for finite momentum states. In cold atom experiments, 
where an external trap is always present, the Fermi sea 
and the homogeneous momentum distribution can be ob- 
served by selectively probing fermions around the centre 
of the trap [si,!!^. 

In this Letter, we have studied the response of a dilute 
two-dimensional Fermi gas with Rashba spin-orbit cou- 
pling to a small overall constant velocity kick. We have 
found that the moving Fermi sea deforms in a non-trivial 
manner due to the non-Galilean nature of the system, 
and is highly degenerate. We have then considered re- 
pulsive interactions at the Hartree-Fock level, and found 
that the ground-state of the system acquires a finite- 
momentum. The Fermi sea becomes deformed beyond 
a critical interaction strength in a continuous fashion, 
which we identified as a possible experimental signature. 
Our results open the path towards the observation of 
finite momentum ground states, constitute the starting 
point for more elaborate treatments of interactions, and 
can be generalized to higher dimensions and more general 
types of spin-orbit coupling. 
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